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Anomalous soft photons in excess of what is expected from electromagnetic bremsstrahlung have 
been observed in association with the production of hadrons, mostly mesons, in high-energy K^p, 
n^p, 7r~p, pp, and e'^e^ collisions. We propose a model for the simultaneous production of anoma- 
lous soft photons and mesons in quantum field theory, in which the meson production arises from 
the oscillation of color charge densities of the quarks of the underlying vacuum in the flux tube. 
As a quark carries both a color charge and an electric charge, the oscillation of the color charge 
densities will be accompanied by the oscillation of electric charge densities, which will in turn lead 
to the simultaneous production of soft photons during the meson production process. How the 
production of these soft photons may explain the anomalous soft photon data will be discussed. 
Further experimental measurements to test the model will be proposed. 

PACS numbers: 24.85.-|-p 25.75.-q 



I. INTRODUCTION 



Anomalous soft photons are soft photons in excess of what is expected from electromagnetic bremsstrahlung. They 
have been observed in conjunction with the production of hadrons, mostly mesons, in K~^p [l|, n^p 0, tt~p 
pp collisions @, and in high-energy e^-e" annihilations in Z° hadronic decay 0-01 • Recent DELPHI measurements 
on the characteristics of the produced hadrons associated with the anomalous soft photon production provide a wealth 
of information on the production process [7l-[l0j. The main features of the anomalous soft photon phenomenon can 
be summarized as follows: 

1. Anomalous soft photons are produced in association with hadron production at high energies. They are absent 
when there is no hadron production [loj . 

2. The anomalous soft photon yield is proportional to the hadron yield. 

3. The transverse momenta of the anomalous soft photons are in the region of many tens of McV/c. 

4. The anomalous soft photon yield increases approximately linearly as the number of neutral or charged produced 
particles increases, but, the yield of anomalous soft photons increases much faster with increasing neutral particle 
multiplicity than with charged particle multiplicity. 

Previously, many different theoretical models have been put forth to explain the anomalous soft photon phenomenon. 
Reviews of the experimental results and theoretical models have been presented [ll|, [l^. T here are models based 
on the assumption of a cold quark gluon plasma [3, boost-invariant classical flux tube [l^, gluon dominance [l5| . 
Unruh-Davies effect ITgI . synchrotron radiation in the stochastic nonperturbative QCD vacuum [l3|, the classical 
string fragmentation [l8|. closed quark-antiquark loop [l^, and ADS/CFT Supcrsymmetric Yang-Mills theory j20j . 
These different proposed models are useful to indicate that to search for the source of the anomalous soft photons, it 
is necessary to go beyond the electromagnetic bremsstrahlung process. While the various models may explain some 
features of the process, the fourth feature listed above from the recent DELPHI observations cannot be explained by 
all existing models [1, @ ; a complete understanding of the basic origin of the anomalous soft photon is still lacking. 

We would like to propose a model for the simultaneous production hadron and anomalous soft photons in the q-q 
string-fragmentation in quantum field theory to explain the anomalous soft photon phenomenon. As described by 



Casher, Kogut, and Suskind [2l|, the production of mesons in such a theory arises from the oscillation of color charge 
densities of the quark vacuum in the flux tube when a quark and an antiquark (or a diquark) pull away from each other 
at high energies. These color charge density oscillations obey the Klein-Gordon equation characterized by the mass 
of the meson [2l| - [3l| . Because a quark carries both a color charge and an electric charge, the underlying dynamical 
motion of the quarks in the vacuum that generate color charge density oscillations will also generate electric charge 
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density oscillations in the flux tube. These color charge density oscillations will lead to the production of photons 
that are clearly additional to those from the electromagnetic bremsstrahlung process. Thus the oscillation of the 
quark densities in the vacuum will lead to the simultaneous color and electric charge density oscillations and will 
subsequently lead to simultaneous and proportional production of mesons and anomalous photons, in agreement with 
the first two features of the anomalous soft photon phenomenon listed in the beginning of this section. 

It is of interest to examine whether the model also leads to results that will be consistent with the remaining features 
of the anomalous soft photon production phenomenon. For such a purpose, we need to know the properties of these 
electric charge density oscillations of the quarks in the flux tube. Will the electric charge density oscillations also obey 
the Klein-Gordon equation characterized by a mass to give rise to stable photons in the flux tube environment? If these 
photons are stable, what are the magnitudes of their masses and how do the masses depend on the quantum numbers 
and other physical properties? If they are produced, how arc they observed experimentally in four-dimensional 
space-time? 

To answer these questions, we start with quarks interacting with both QCD4 and QED4 interactions in four- 
dimensional space-time in Section II. We specialize to the case of the flux tube formation for high energy particle 
production processes under longitudinal dominance and transverse confinement. The system can then be approx- 
imately compactified into a QCD2xQED2 system in two-dimensional space-time, with the coupling constants in 
different space-time dimensions related by the flux tube radius. In Section, III, we examine the non-Abclian bosoniza- 
tion of the QCD2xQED2 system. We find stable QCD2 and QED2 bosons for quarks with two flavors arising from 
the density oscillations of the quarks in the flux tube. These bosons can be identified as QCD2 mesons and QED2 
photons. The boson masses are then expressed as a function of the coupling constants and the quark transverse 
mass. In Section IV, we estimate the coupling constants and the boson masses for the case of the hadronic decay. 
As QCD2 mesons and QED2 photons are stable in the flux tube environment, we can infer from the quantum field 
theory description of the particle production process in Ref. [2l| that these QCD2 mesons and QED2 photons will be 
produced simultaneously in q-q string fragmentation. The present model for the production of anomalous soft photons 
can also be conveniently called the QED2 photon model. In Section V, we discuss the adiabatic decompactification 
of produced photons and mesons from two-dimensional space-time to particles in four-dimensional space-time. In 
Section VI, we investigate how the QED2 model of photon production may explain the experimental anomalous soft 
photon transverse momentum distributions. In Section VII, we examine the rates of QED2 meson and QED2 photon 
production and the correlation of the soft photon yield with charge and neutral particle multiplicities. In Section 
VIII, we suggest future experimental measurements to test the QED2 photon model. In Section IX, we present our 
conclusions and discussions. 



II. FLUX TUBE ENVIRONMENT IN HIGH ENEGY PARTICLE PRODUCTION PROCESSES 



We wish to investigate the process of soft photon production in association with hadron production, in which the 
produced hadrons consist mostly of mesons. In the quantum field theory description of meson production process as 
analogous to the particle production process in quantum electrodynamics in two-dimensions (QED2), mesons that 
are stable within the theory will beproduced along the string, when a quark and an antiquark at the two ends of the 
string pull apart at high energies |2ll425l | . The rapidity distribution of these produced mesons exhibits the property 
of boost invariance in the limit of infinite energies pTI . \2l\ . For a finite energy system, the boost-invariant solution 
turns naturally into a rapidity plateau, whose width increases with energy as \n{^/s) [281 [29|. 

The q-q string is an idealization of a flux tube with a transverse profile, which reveals itself as the transverse 
momentum distribution of the produced particles [32| . Exp eriment ally , the presence of a flux tube is evidenced by the 
limiting average transverse momentum and a rapidity plateau [isl . uu . [26l - [3^ as in high-energy e+-e~ annihilations 
[H-ll^l and pp coflisions HI]. 

To investigate the simultaneous production of QCD and QED quanta in the flux tube fragmentation, we study 
quarks interacting with both QCD and QED interactions, in circumstances leading to the formation of the flux tube. 
It becomes convenient to consider the U(3) group which breaks up into the color SU(3) and the electromagnetic 
U(l) subgroups. The SU(3) and U(l) subgroups differ in their coupling constants and communicative properties. We 
introduce the generator t° for the U(l) subgroup, 

1 /I o\ 

i° - 1 , (1) 

which adds on to the eight generators of the SU(3) subgroup, {t^, ■■■t^}, to form the nine generators of the U(3) group. 
They satisfy tr{t"i^} = (5"^/2 for a, /3 = 0, 1, .., 8. 
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Limiting our consideration to quarks with two light flavors, we examine the QCD4xQED4 system in four- 
dimensional space-time x'', with /x=0,l,2,3. The dynamical variables are the quark fields, ip'j, and the U(3) gauge 
fields, Ai,=A'j^t°', where a is the color index with a=l,2,3, / is the fiavor index with f=u, d, and a is the U(3) generator 
index with q;=0,1,..,8. The coupling constants depend on a and / and are given explicitly by 

gi' -'^ = <?1'' "'^ = 5QCD4, for QCD, (2) 

= -Cm = -Q«eQED4, 5° = -e<i = -Q<ieQED4 for QED, (3) 

with Qu = 2/3, and Qd = —1/3. We use the convention of summation over repeated indices, but the summation 
symbol and indices are occasionally written out explicitly to avoid ambiguities. For brevity of notations, the indices 
a, /, and a in various quantities are implicitly understood except when they are needed. For example, the 6-color 
component of (g^^V)'' written explicitly is Ea=o,...,8 Ea=i,2,3 

The transverse confinement of the flux tube can be represented by quarks moving in a transverse scalar field m(r) 
where m{r) = S'(r)-|-(currcnt quark mass m^) and S{r) is the confining scalar interaction arising from nonperturbative 
QCD. The equation of motion of the quark field ip is 

{t]p-m{r)}i'^0, (4) 

where 

iP^^^n^^^^ip^+gA,). (5) 
The equation of motion for the gauge field A^^ is 

D.F^'' = d^F^-' - ig[A^, F^^] = gf, (6) 

where 

F^u = df,A^ - d^A^ - ig[A^„ A,,] , (7) 

F,.^F^X: f^ffX: (8) 

j7«=2trV5/7''t«V/- (9) 

Because of the commutative properties of the generator, the commutator terms in Eqs. ([B]) and ([7]) give zero 
contributions for QED. This set of equations describe the coupling of the QCD and QED gauge fields and the 
quark fields, including those in the quark vacuum. The self-consistent coupling of these fields lead to a problem of 
great complexity. Fortunately, they can be simplified under the conditions of longitudinal dominance and transverse 
confinement that exist in the particle production environment at high energies [3l|. 

This set of coupled equations in four dimensional space-time for QCD4xQED4 with the U(3) interaction has the 
same mathematical structure as our previous set for QCD4 with the SU(3) interaction [3l[. As the approximate com- 
pactification depends on the separation of the transverse and the longitudinal degrees of freedom and is independent 
of the nature of the underlying gauge group, there should be similar approximate compactification of QED4xQED4 
into QCD2xQED2. Wc shall briefly summarize the salient points leading to such an approximate compactification 

In the problem of particle production at high energies leading to the formation and fragmentation of the fiux tube, 
we can focus our attention on the self-consistent coupling of the quarks and the gauge fields A^ inside the tube. 
In these high-energy processes, the momentum scales for longitudinal dynamical motion of the leading q and q as 
well as those of quarks in the underlying vacuum are much greater than the momentum scales for their transverse 
motion such that \v'^\ ^ jw^j, jw^j, where v is a typical quark velocity. In the Lorentz gauge, the associated gauge 
field is proportional to (l,v). Under the dominance of the longitudinal motion over the transverse motion in 
string fragmentation, \A^\, \ A^\ ^ |^^|, \ A?'\. Hence, inside the flux tube A^ and A? can be approximately neglected 
in comparison with the magnitudes of A^ or A^ . It is further reasonable to assume that the gauge fields A^ and A'^ 
in the interior of the tube depend only weakly on the transverse coordinates r = [x^^x^). It is then meaningful to 
investigate these fields inside the tube by averaging them over the transverse profile of the fiux tube. After such an 



averaging, and inside the tube can be considered as a function of {x^ , x^) only. As a consequence, the equation 
of motion 1^ for the quarks becomes 

{7°no + 7 + 7^P2 + 7^n3 - m(r) } V = 0. (10) 

We write the quark field ip{x) in terms of the longitudinal fields f± and transverse fields Gi^2 with spinors as [l^ 

i;{x) = [Gi(r)6 - G2(rK2]/+(x",x3) + [Gi(r)6 + G2{r)Uf^{x'' , x'), (11) 

where 





= h h 6 = n h €3 = I _i I , and ^4 = q | • (12) 

Vl. 

Working out the Dirac matrices in we obtain the following set of coupled equations of motion 

[(n"-n3)(n° + n3)-m2(r)-p|]Gi/_ = -[{p' -2p^)m{r)]G2f^, (13a) 

[(n°-n3)(nO + n3)-m2(r)-p|]G2/- = [(/ +zp2)m(r)]Gi/_, (13b) 

[(nO + n3)(n"-n3)-m2(r)-p|]Gi/+ = -[{p' ~ zp^)m{r)]G2f+, (13c) 

[(n" + n3)(nO-n3)-m2(r)-p|]G2/+ = [(/ +zp2)m(r)]Gi/+. (13d) 

By the method of the separation of variables, we introduce the eigenvalue (the quark transverse mass) uit for transverse 
motion, 

[p| + m2(r) -m|]Gi(r) = [{p^ - ip'')m{r)]G2{r), (14a) 

[p| + m2(r) -m|]G2(r) = -[(p^ + ip2)TO(r)]Gi(r), (14b) 

and obtain the coupled equations for longitudinal motion, 

[(n0-n3)(n" + n3)-m|]/_(x°,x3) = o, (i5a) 

[(n° + n3)(n"-n3)-774]/+(x°,x3) = o. (i5b) 



If we introduce the two-dimensional Dirac spinor -020 as 



^2D =(/_): (16) 



and the 2-dimensional gamma matrices as [4l|, |48 



72D ^ (^1 Oj ' =[^1 j ' '^2D72D = 72D = (^Q -1 j ' (^'^^ 

then Eqs. (|15ap and (|15b[) can be rewritten as the Dirac equation 

{72Dbo + 9A0) + 7|d(p3 + g^s) - mr} ^^20 = 0, (18) 

which is the equation of motion for a quark in two-dimensional gauge fields of Aq and A3, except that the coupling 
constants g are those in four-dimensional space-time, g^D, and the gauge fields v4o(4d) and ^3(40) are those determined 
from a four-dimensional current source j^Y) given by Eqs. ^ and © involving four-dimensional quark fields ipAn- It is 
necessary to renormalize the coupling constants and use quantities determined from two-dimensional source currents 
and fields, to bring it to the proper two-dimensional space-time form. 

Utilizing the result of Eq. ([T51) for the quark field and using the quark wave function ([TI]) , the set of equations of 
motion Q-® along the longitudinal direction can be cast into the forms of quarks and gauge fields interacting with 
the QCD and QED in two-dimensional space time of (a;°, x^), by transversely averaging Eq. ([5]) over the profile of the 
flux tube and by relating the coupling constants with the renormalization |3l| 

3lD=94D((|Gi(r)p + |G2(r)n)^. (19) 



Such a renormalization yields a two-dimensional coupling constant 320 that possesses the dimension of a mass. 
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After the coupling constant rcnormalization, the equations of motion for the quark and the gauge fields in the 
longitudinal and time directions are as given in the same form as Eqs. ([I])-® in two-dimensional space-time a;''(/i = 
0, 3) with ^/'4D replaced by "020 , 74d ganima matrices replaced by the 2-dimensional gamma matrices 72Dj the quantity 
m(r) replaced by the quark transverse mass tot, and gauge fields A^(4d) limited to /x = 0,3 and replaced by Af^(^2B) 
determined from the Maxwell Equation D^F^^ = 52DJ2D' with two-dimensional currents j'jd that arise from %I]2D- 

We can get an approximate relation between 920 and 940 by considering the case of a uniform transverse flux tube 
profile with a transverse radius i?T, 

(|Gi(r)p + |G2(r)p) ^ e(i?T - |r|)/7ri?|,. (20) 

The coupling constants in two-dimensional space-time and 4-dimensional space-time are then related approximately 
by mi 

Such a result is expected from dimensional analysis, where the length scale in going from a tube to a string involves 
only the flux tube radius. The above relationship between the coupling constants in different space-time dimensions 
will be used later to estimate the boson masses. 



III. BOSONIZATION OF QCD2xQED2 FOR QUARKS WITH TWO FLAVORS 

Under the longitudinal dominance and transverse confinement, the QCD4xQED4 system can be approximate 
compactified as the QCD2xQED2 system with a quark transverse mass toj-. The flux tube becomes the arena for 
the quarks in the underlying vacuum to interact self-consistently with the QCD and QED gauge flelds. We shall 
henceforth work with QCD2xQED2 in two-dimensional space-time. For brevity of notation, the two-dimensional 
designation of various quantities will be understood in what follows. The Lagrangian density for QCD2xQED2 that 
corresponds to the two-dimensional version of Eqs. (jl])-® is 

£ - V5[7^(^9^ + gA^) ~ mT]iJ - ^F^.F'^". (22) 

As in the Section II, the color index a, the flavor index /, and the U(3) generator index a are implicitly understood, 
and the summation convention is used. 

We wish to search for bound states arising from the density oscillations of the color and electric charges of the 
quarks in QCD2xQED2 in the strong coupling limit, in which the strength of the QCD2 interaction is much greater 
than the quark mass. The best method to search for bound states is by bosonization in which bosons are bound and 
nearly free, with residual sine-Gordon interactions that depend on the quark mass [2^.[39|-[58j. 

The U(3) gauge interactions under consideration contains the non-Abelian color SU(3) interactions. Consequently 
the bosonization of the color degrees of freedom should be carried out according to the method of non-Abelian 
bosonization which preserves the gauge group symmetry (4l| . 

While we use non-Abelian bosonization for the U(3) interactions, we shall use the Abelian bosonization for the 
flavor degrees of freedom. This involves keeping the flavor labels in the bosonization without using the flavor group 
symmetry. Although the Abelian bosonization in the flavor sector obscures the isospin symmetry in QCD, the QCD 
isospin symmetry is still present. It can be recovered by complicated non- linear general isospin transformations 

[Hiol. 

As in any method of bosonization, the non-Abelian method will succeed for systems that contain stable and 
bound boson states with relatively weak residual interactions. Thus, not all the degrees of freedom available to the 
bosonization technique will lead to good boson states with these desirable properties. For example, some of the 
bosonization degrees of freedom in color SU(3) may correspond to bosonic excitations into colored objects of two- 
fermion complexes and may not give rise to stable bosons. It is important to judiciously search for those boson degrees 
of freedom that will eventually lead to stable and bound bosons. 

Keeping this perspective in our mind, we can examine the non-Abelian bosonization of the system under the U(3) 
interactions. The non-Abelian bosonization program consists of introducing boson flelds to describe an element u of 
the U(3) group and showing subsequently that these boson flelds lead to stable bosons with finite or zero masses. 

In the non-Abelian bosonization, the current j± in the light-cone coordinates, x^ = {x^ ± x^)/\/2, is bosonized as 

m 

j+ = {t/27T)u~\d+u), (23a) 
j_ = ~{i/2Tr){d_u)u-\ (23b) 
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An clement of the U(l) subgroup of the U(3) group ean be represented by the boson field (fp 



u 



exp{i2V^</>°i°}. (24) 



Such a bosonization poses no problem as it is an Abelian subgroup. It will lead to a stable boson as in Schwinger's 
QED2. 

To carry out the bosonization of the color SU(3) subgroup, we need to introduce boson fields to describe an element 
u of SU(3). There are eight generators which provides eight degrees of freedom. We may naively think that for 
the non-Abelian bosonization of SU(3), we should introduce eight boson fields to describe u by 



7i = exp{^2^AF^ (25) 

Q = l 

However, a general variation of the element 6u/Sx^ will lead to quantities that in general do not commute with u 
and u^^, resulting in j± currents in Eqs. (j23p that are complicated non-linear admixtures of the boson fields 0". It 
will be difficult to look for stable boson states with these currents. 

We can guide us to a situation that has a greater chance of finding stable bosons by examining the bosonization 
problem from a different viewpoint. We can pick a unit generator n = {n^, n^, .., n^} oriented in any direction of the 
eight-dimensional a-space and can describe an SU(3) group element u by an amplitude and the unit vector n, 

8 

u cxp{^2^/¥<^ n"i"}. (26) 

The boson field (j) describes one degree of freedom, and the direction cosines {n",a = 1,..,8} of the unit vector 
n describe the other seven degrees of freedom. A variation of the amplitude (f> in u while keeping the unit vector 
orientation fixed will lead to a variation of 5u/Sx^ that will commute with u and u^^ in the bosonization formula 
(1231), as in the case with an Abelian group element. It will lead to simple currents and stable QCD2 bosons with well 
defined masses, which will need to be consistent with experimental QCD meson data. On the other hand, a variation 
of 6u/6x^ in any of the other seven orientation angles of the unit vector n will lead to Su/Sx^ quantities along other 
t" directions. These variations of Su/6x^ will not in general commute with u or u~^. They will lead to j± currents 
that are complicated non-linear functions of the eight degrees of freedom. We are therefore well advised to search 
for stable bosons by varying only the amplitude of the (j) field, keeping the orientation of the unit vector fixed, and 
forgoing the other seven orientation degrees of freedom. 

As a unit vector n in any orientation can be rotated to the first axis along the direction by an orthogonal 
transformation in the a-space, we can consider the unit vector n to lie along the direction without a loss of 
generality. For the U(3) group. The appropriate bosonization program that will eventually lead to stable bosons is to 
limit the consideration to only the and (j>^ degrees of freedom. We are therefore justified to bosonize an element u 
of the U(3) group as 



exp{i2V^^(^"r}. (27) 



a=0 



From Eqs. (P^a)) and ([23b)) . we obtain then 

1 

/TT 



1 1 



a=0 



where we have written out the flavor index explicitly. The gauge fields can be easily obtained by using the = 
gauge for which terms involving the commutators in Eqs. ([6]) and (jT)) vanish. The Maxwell equation becomes 



-dlA+ = -gj+ (29) 

and the solution is 



= + . (30) 



The interaction energy becomes 



H, = \l dx^jlpl = ^ /^-- E( E 5/0/)^- (31) 



=0 f=u,< 
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The kinetic energy term of the Lagrangian density, ij}^^ id^ijj , can be bosonized as [41 1 



f=u,d 



as the Wcss-Zumino term for u in the form of Eq. ([27|) gives no contribution. Eq. ((27)) then leads to 

f=u,d a=0 

The mass term involves the scalar density ip4' which can be bosonized as 



(32) 



(33) 



2^ 



271 



Uf + Uj:' 



f=u,d 



tiN^ E tr[cos(2V^E'^/^")]' 

f=u,d a=Q 



(34) 



where 7 = 0.5772 is the Euler constant, N^^ is normal ordering with respect to the mass scale n for the problem in 
question. It is easy to show that 



tr[cos(2^/^E = 2 cos(2V^0/) cos(2yV40/)- 



(35) 



We shall not examine the zero mode and the theta vacuum in the present exploratory study. We obtain the Hamiltonian 
density 



a—0 ^ f—u,a 



f=uA 



27r 



2N^, E cos(2V^0/)cos(2yV4'/>/)- 



(36) 



In the flavor sector, the up quark has isospin quantum numbers (/, /3)=(l/2, 1/2) and the down quark has 
(/, /3)=(l/2, — 1/2). The up and down quarks combine to form the isoscalar (/, /3) = (0, 0) state and the isovector 
1=1 states, which split into three components with /3 = (1,0,— 1). Because the quark electric charge Qf depends on 
the flavor quantum number, there is no isospin symmetry for QED2, and the four states split apart. We shall focus 
our attention only on the isoscalar (/,/3)=(0,0) QED2 state and the isovector (/,/3)=(l,0) QED2 state. The other 
two (/, /3) = (1, ±1) QED2 states involve composite constituents with like electric charges and are unlikely to be stable 
in the electromagnetic sector. For brevity of nomenclature, we shall refer to the isovector (/, /3)=(1, 0) photon simply 
as isovector photon, with the qualifying specification '/3=0' implicitly understood. 

In QCD with two flavors, the isospin symmetry remains a good symmetry, which is weakly broken by a small 
difference between the up and down quark masses. Thus, the QCD quark-antiquark meson states are specified by 
isospin quantum numbers / with nearly degenerate 2/+ 1 members of different /a components. The knowledge of the 
location of the (/, /3)=(1,0) QCD state allows us to infer the locations of the other two QCD (/, J3)=(l, ±1) states. 

We can construct the fields for the isospin (/, = 0) states, for up and down quark fields moving in phase or 
out of phase, 



We can also construct the corresponding isospin canonical momenta 

= i= [uz + (-i)^ns] 



(37) 



(38) 



The Hamiltonian density in terms of boson fields of different isospin quantum numbers / and the same /3 = is 

1 1 



EE 



J(n?)^ 



(39) 
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where V{{(t)"}) — Vint{{4'f}) + Vm({0/ }) with the interaction energy 

-(-1 
^/27^ 



VUm) = i;\£^-^^^^ro] , (40) 



\I=0 



and the quark mass term 



27r 



(41) 



We can get the gross features of the system by expanding the potential about the minimum located at (/)g = and 
0" = 0. Evaluating the second derivatives of the potential at the potential minimum, we obtain the mass square 
(A'/"(.2Q-))^ of stable boson quanta for q;=0,1, 



(m;(2D)) 



v{{m 



9Z + 



27r / 3 - a 



e'^TOTM- (42) 



The Hamiltonian density ([39|) represents a QCD2 and QED2 system of isoscalar and isovector boson fields 0" whose 
field quanta acquire the mass Af"^2D)' '^here a = for QED2 and a = 1 for QCD2. As the boson field 0" is related 
to the gauge field through Eqs. (29) and (31), the quanta of are also the quanta of the gauge fields A+. The 
QCD2 bosons and QED2 bosons can be appropriately called QCD2 mesons and QED2 photons respectively. 

Because the righthand side of Eq. (|42|) is a non-negative quantity with {Mf^^j^^Y ^ 0> these QCD2 mesons and 
QED2 photons are stable bosons. They acquires a mass because a gauge field oscillation leads to a quark density 
oscillation, and through the Maxwell equation the quark density oscillation in turn leads to a gauge field oscillation, 
which in turn modifies the quark density oscillation. The self-consistency of gauge field oscillations and the induced 
quark density oscillations lead to an equation of motion for the gauge field oscillation in the form of a Klein-Gordon 
equation characterized by a mass. 

Our result of the boson masses in Eq. ()42|) represents a QCD2xQED2 generalization of previous results in p2l|.(40l|- 
[5^, where QED2 and QCD2 have been examined separately. In the massless quark limit, the QCD2 and QED2 
boson masses are given by jg" -I- (— l)^5^|/-\/27r. In this limit, the QCD2 masses are the same as what one obtains by 
using QED2 and replacing the electric charges in QED2 with the color charges in QCD as in [22j. This equivalence 
of the Abelian QED2 solution and the non-Abelian QCD2 solution in the massless limit arises because our judicious 
search for stable QCD2 bosons in the non-Abelian bosonization of SU(3) requires the variation of only the amplitude 
(f) while the orientation of n in Eq. ([25]) is held fixed. The non-Abelian bosonization in QCD2 that results in stable 



QCD mesons is in effect Abelian in nature. This explains why previous Abelian QED2 results of boson masses [22 1 
and string fragmentation [2l| can be applied to the non-Abelian QCD problems by replacing the electric charges in 
QED2 with the color charges in QCD. 

Previously, Abelian-type solutions were obtained for multiflavor QCD2 mesons using non-Abelian bosonization for 



both the color and flavor degrees of freedom in the large Nf limit 52|. The mass of the single massive boson in the 

"60I ]. Our QCD2xQED2 analysis here indicates 



massless quark limit was found to be M^^qqy^2) ~9qcd2 \/ Z"" [Si 
that Abelian-type solutions exist also for QCD2 mesons for quarks with two flavors, and is not limited to the large Nf 
limit. Our mass of the QCD2 isoscalar meson in the massless quark limit is -^^Q(QCD2)~ffQcr>2 ■v/^/tt, which matches 
the mass of the massive boson of [s^ - lssl . l60j for Nf=2. Thus, by using the non-Abelian bosonization in QCD2 but 
Abelian bosonization in the flavor degrees of freedom in the present treatment, the solutions of [52l - [55| in the large 
Nf limit can be extended down to iV/=2. 

In the massless quark limit (for mTi^=0 in this case), the QCD2xQED2 bosons are free. With a finite value of 
rriTfJ-, they interact with a sine-Gordon residual interaction whose strength depends on rriT^- The present treatment 
places the QED2 mesons and the QCD2 photons on a parallel footing and allows the mutual interaction between 
QCD2 mesons and QED2 photons. To exhibit the mutual interaction, it is instructive to expand the quark mass term 
in powers of 0". Up to the fourth order in (/)", we obtain 

VMf)) = \Y.jZamf + \jZ E + jEEczz'(0?f(0}f (43) 

Q=0 7=0 Q=0/./'=0 Q=0 7=0 
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where 

,44a) 

b7r = --r^e'<mT^^, (44b) 
(3 - ay 

7T 

Cjj' = ——eJjnTfJ^. (44c) 
o 

Here the af coefficients give the contribution to (M^i^^d))'^ from the quark mass term in Eq. (|^^ . The bfj, coefficients 
give the interaction between bosons of the same type a, and c//' give the interaction between QCD2 mesons and 
QED2 photons. The negative signs of the b and c coefficients indicate that the interaction between the bosons are 
attractive in nature. 

Previously. Coleman obtained the correction to the QED2 boson masses arising from a non-zero quark mass, using 
the method of renormal-ordering. The mass correction was also obtained by examining QED2 on a circle [sij , near- 
light-cone coordinates [i^ , and mass-perturbation theory [s^ . We shall not consider these refinements and contend 
ourselves with the estimate using the second derivatives of the potential V{{(f>f}) in the present exploratory study. 



IV. QCD2 MESON AND QED2 PHOTON MASSES FOR QUARKS WITH TWO FLAVORS 

We consider first the boson masses in the massless quark limit because they represent well-defined references. In 
this limit, the boson masses depend only on the coupling constants which acquire the dimension of a mass as a result 
of the compactification. They depend on the flux tube radius as given by Eq. (PT|) [30l . Isi} . For QCD in the flux tube, 
the QCD2 coupling constant is given by 

2 2 

2 gQCD4 _ gQCD4 4 _ Aa^ , , 

where as = 5qcd4/4''' is the strong interaction coupling constant. Similarly, for QED2 in the flux tube, the QED2 
coupling constant is given by 

2 _ ^QED4 _ gQED4 4 _ 4a 

where a = eQ^j^^/Air = 1/137 is the fine structure constant. The magnitude of the flux tube radius Rt is revealed 
by the root-mean-squared transverse momentum of produced hadrons (mostly pions) as 

Rt , (47) 

which empirically is slightly energy-dependent [3l| . We shall focus our attention on the case of particle production 
in high energ y e+-e~ annihilations in the hadronic decay of Z^. The measurement of the 7r° spectra in Z° hadronic 
decay gives \/ {pt)tt = 0.56 GeV in the reaction plane [61| and thus the flux tube has a radius Rt ^ 0.35 fm. For 
the strong cou plin g constant at this energy, we shall take a^ = 0.316, which leads from Eq. to the string tension 
coefRcient [30l.l3l| 

& = gQCD2/2-0.2GcV2, (48) 

and 

ffQCD2 = 0.632 GeV. (49) 

From Eq. (|46p . the QED2 electromagnetic coupling constant has the value 

eQED2 - 0.096 GeV. (50) 

With these coupling constants (3i=3d=5QCD2, ff2=-Q«eQED2, and 5°=-QdeQED2), the values of QCD2 and QED2 
boson masses in the massless quark limit are shown in Table I. One observes that QCD2 for quarks with two flavors 
gives a massless pion in the massless quark limit, in agreement with the concept of the pion being a Goldstone boson 
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TABLE I: QED2 and QCD2 boson masses obtained with _Rt=0.35 fm and gQCD2=2&=0.4 GeV^. 





QCD2 


QED2 


Coupling Constant 


(7QCD2=632.5 MeV 


eQED2=96 MeV 


massless quarks 


isoscalar boson mass Mo(2d) 
isovector boson mass Mi(2d) 


504.6 MeV 



12.8 MeV 
38.4 MeV 


mT=400 MeV 


isoscalar boson mass Mo(2d) 
isovector boson mass Mi(2d) 


734.6 MeV 
533.8 MeV 




mT= 400 MeV 
H=mq=0{l MeV) 


isoscalar boson mass Mo(2d) 
isovector boson mass Mi(2d) 




0(25.3 MeV) 
0(44.1 MeV) 



in the standard QCD theory. The isovector QCD2 meson lies lower than the isoscalar QCD2 meson at 504 MeV, 
whereas the ordering is opposite for the QED2 photons, with an isoscalar QED2 photon at 12.8 MeV and an isovector 
QED2 photon at 38.4 MeV. These QED2 photons lie in the region of observed anomalous soft photons. 

Equation (|42p indicates that the boson masses depend on four mass scales: 5qcd2, eQED2, tot, and /i. In addition 
to the coupling constants we have just discussed, we need to specify the values of the transverse mass mx and the mass 
scale fi. The discussions in the Section II indicate that as quarks resides in the flux tube environment, they acquire 
a transverse mass niT- The presence of the factor hit in Eq. (j42p takes into account the effects of non-perturbative 
chiral symmetry breaking and transverse confinement that lead to the formation of the flux tube. Because a pion is 
a quark- antiquark composite, we can estimate the quark transverse mass niT from the pion transverse momentum, 
rriT ^ \J {pt)-k/'^- For Z° hadronic decay, {px)-K = 0.56 GeV and we have rriT ~ 0.4 GeV. 

The boson masses depend also on the mass scale /x, which arises from the bosonization of the scalar density ^ij} 
as given in Eq. p4|) . The scalar density ij^i}) diverges in perturbation theory and has to be renormalized such that 
(^'ipip'^=0 in a free theory. It will need to be rcnormal-ordered again in an interacting theory [2^ . The scalar density 
and the corresponding mass scale therefore depends on the interaction. The dependence of the scalar density on the 
interaction is also evidenced by the fact that the scalar density "ipip in Eq. (|34p can be expanded in terms of i° and t^, 
each of which is the generator of a different interaction associated with a different coupling constant. For the strong 
interaction of QCD, confinement and chiral symmetry breaking dominate and lead to a transverse mass rriT that is 
much greater than the current quark mass. It is reasonable to take the mass scale fi in QCD to be the same as the 
quark transverse mass rriT characterizing the flux tube transverse confinement and the presence of chiral symmetry 
breaking. Meson masses in QCD2 calculated with the mass scale ^ = mr = 0.4 GeV is given in Table I. It gives a 
QCD2 isovector meson mass of 0.534 GeV and a QCD2 isoscalar meson mass of 0.735 GeV in the flux tube. 

For a theory with a relatively weak interaction such as QED, the scalar density 'ipip that diverges in perturbation 
theory has to be renormalized in a nearly-free field in which the quark energy is just the current quark masses. The 
mass scale /i for QED2 should therefore be the QED current quark mass niq appropriate for a nearly-free theory. 
The current quark mass ruq associated with perturbative QCD has the value of 1.5-6 MeV [62|. The current quark 
mass associated with perturbative QED is not known and presumably is of the same order of an MeV. For lack of a 
more definitive determination, we shall take /i ~ 1 MeV to calculate the orders of magnitude of the QED2 photon 
masses. The values of the QED2 boson masses obtained with n = 1 MeV are given in Table I, which gives an isoscalar 
photon of order 25 MeV and an isovector photon of order 44 MeV. They fall within the same order of magnitude of 
the transverse momenta of anomalous soft photons. 

V. ADIABATIC DECOMPACTIFICATION OF BOSONS FROM TWO-DIMENSIONAL TO 

FOUR-DIMENSIONAL SPACE-TIME 

We have thus found that in the system of quarks with two flavors, the boson quanta of QCD2 and QED2 are 
stable with masses that depend on the isospin quan tum numbers. We can therefore infer from the quantum field 
theory description of particle production in Ref. [2l| that these QCD2 mesons and QED2 photons will be produced 
simultaneously in the same process of q-q string fragmentation, when a quark pulls away from an interacting antiquark 
at high energies. 

After a particle is produced in two-dimensional space-time how docs it dccompactify in the four-dimensional space- 
time? An appropriate way to describe the decompactification is to identify the mass of the particle in the two- 
dimensional theory as the transverse mass of the particle in four-dimensional space-time. This clearly works in the 
case of the quark. In reverting back into the four-dimensional space-time, the 'quark mass' of mx in two-dimensional 
space-time reverts back into the transverse mass of the quark in four-dimensional space-time. 



11 



After a boson is produced and the system expands longitudinally, the interaction between the produced bosons 
weakens. The produced boson will subsequently emerge from the production region out to the non-interacting region 
and will obey the mass shell condition. We can consider a produced boson of mass M"^2d) isospin / and type a 
in two-dimensional space-time. The kinematic variables of the boson are i?(2D) Pz(2D): which obey the mass shell 
condition E'^2D) ~P'z {2D) ~^ (-^ T{2D))'^ ■ ^^'^ four-dimensional space-time, kinematic variables of this particle are £'(4D)j 
Pz(4D)i and Pt(4D)- The kinematic variables satisfy the mass shell condition £'^\p)=p^(4j3^+Py(4j-,)-|-(M"^4j-,^)2, where 
^T{4D) '^^ ^^^^ mass of the particle in four-dimensional space-time. 

In the emergence of the boson from two-dimensional space time to four-dimensional space time, we envisage an adi- 
abatic transverse expansion from the two-dimensional flux tube to four-dimensional space-time. The adiabatic trans- 
verse expansion involves no change of the particle energy and particle longitudinal momentum so that i?(2D)=^(4D) s-nd 
Pz(2D)=Pz{4D)- Consequently, the mass shell conditions of the boson give (-M7(2D))^=Pt(4D) + (-^/(4D))^=(-^/t)^; '^ith 
the boson mass Mf^2T:,) two-dimensional space-time turning into the boson transverse mass in four-dimensional 
space-time. 

We can test the consistency of such a correspondence for the production of mesons. We consider the production of 
an isovector meson, which is a pion. Experimentally, a pion is produced with an average \/Jpt)^ ^ 0.56 GeV 
the hadronic decay. Thus, the experimental (average) isovector meson (pion) transverse mass is 

M^V = V (0-14)2 ^ (0.56)2 ^ 0.579 GeV. (51) 

We consider next the production of an isoscalar meson, which is the rj meson with a rest mass A/^ = 0.547 GeV. 
Experimentally, the observed average transverse momentum of a meson increases with the meson mass. The average 
transverse momentum of rj has however not been measured. As the rj meson has approximately the same mass as 
a kaon whose average transverse momentum has been measured, the a verage t ransverse momentum of the r] meson 
should be of the order of the kaon average transverse momentum of \/ (p'^)k ^ 0.616 GeV [ssf. Upon taking this 
estimate to be the isoscalar meson average transverse momentum, the experimental (extrapolated) isoscalar meson 
(77 meson) average transverse mass is 

A/q'V = V(0-547)2 + (0.616)2 GeV - 0.824 GeV. (52) 

We can compare the experimental average transverse masses of isovector and isoscalar mesons in Eqs. ([ST]) and ((52|) 
with the theoretical QCD2 meson masses in two-dimensional space-time in Table I, which gives Af7r(QCD2) =0.534 GeV, 
and Afjj(QC]32) =0.735 GeV. We find that there is approximate agreement of the experimental meson average transverse 
masses in four-dimensional space-time with the QCD2 meson masses in two-dimensional theory, within about 10-15%. 
This approximate agreement lends support to the identification of the mass of a stable boson in QCD2xQED2 as the 
(average) transverse mass of the boson in four-dimensional space-time. 

In the case of QED photons, the rest mass of the photon is zero in four-dimensional space time. Hence, the QED2 
photon mass in two-dimensional space-time can be identified as the photon transverse momentum in four-dimensional 
space-time. 



61| for 



VI. ANOMALOUS SOFT PHOTON TRANSVERSE MOMENTUM DISTRIBUTION 



We shall explore how the model of simultaneous meson and photon production in the string fragmentation process 
may explain the anomalous soft photon phenomenon. In e^-e^ annihilations or hadron-hadron collisions, q-q strings 
or q'-(diquark) strings will be formed with a quark and an antiquark (or a diquark) pulling apart at the two ends 
of each string. As QCD2 mesons and QED2 photons are found to be stable bosons in QCD2x QED2, we can infer 
from the quantum field theory of particle production as described by Casher, Kogut, and Suskind 2l| that QCD2 
mesons and QED2 photons will be produced simultaneously in the same process of q-q string fragmentation, when the 
quark pulls away from the antiquark (or diquark) at high energies. The simultaneous production from the same string 
explains why anomalous soft photons are present only in association with hadron production and why the number 
of produced mesons and QED2 photons are proportional to each other, in agreement in the first two features of the 
anomalous soft photon phenomenon listed in the Introduction. We shall examine whether the QED2 photon model 
can explain the transverse momentum distribution in this section and the correlation of the soft photon yield with 
hadron production properties in the next section. 

According to the Schwinger's mechanism [63l - l66j . the probability of particle production is an exponential function 
of the square of the transverse mass of the produced particle. For massless photons, the photon transverse mass 
is the photon transverse momentum. It is therefore reasonable to assume the transverse momentum distribution 
of each photon component to be a Gaussian with an average root-mean-squared transverse momentum Mjj, given 
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FIG. 1: (Color Online) Comparison of normalized anomalous soft photon dN'' /dpr per pp collision at 450 GeV/c after 
subtracting the bremsstrahlung contributions [gJ with calculated distributions based on three components of anomalous 
soft photons at 17, 45, and 5 MeV. The solid curve is the total theoretical dN'' /dpr distribution, and the other curves are the 
separate contributions from the three different components. 



by the QED2 photon mass. In the measurement of the soft photon transverse momentum distribution in e+-e^ 
annihilations, the determination of the orientation of jet axis has an uncertainty of ~50 mrad, corresponding to 
a root- mean-square uncertainty of ApT —E^lS.6^\Q MeV for the small pt region Q. In the measurement of the 
transverse momentum distribution in pp collisions, the uncertainty in angular measurements is AS'-^^IO mrad, w^hich 
corresponds to lS.pT^2 MeV for the small pr region Q. These uncertainties in the determination of the angles need 
to be folded into theoretical calculations in order to compare with experimental data. If one assumes a Gaussian 
distribution of the transverse coordinates in the angular determination, the folding of two Gaussian distributions 
leads to a Gaussian distribution with a standard deviation square of [(Mjj,)^ -I- (Apt)^]/2. 

Based on our theoretical results in Table I, we expect that there will be two components in the soft photon transverse 
momentum spectrum. We therefore parametrize the experimental anomalous soft photon transverse momentum 
distribution as the sum of the two normalized Gaussian components of isoscalar and isovector photons, each of which 
has an average root-mean-squared transverse momentum given by {M]rpY + {^PtY ■ We search for distributions 
characterized by a transverse mass in the region of M^rp^25 MeV for the isoscalar photon and around M^rp^^iA McV 
for the isovector photon component. However, upon a careful examination of the transverse momentum distribution 
of the anomalous soft photons in pp collisions [6|, we find that in addition to these isoscalar and isovector photon 
components, the pp data appears to contain an additional lower momentum component characterized by a small mass 
Mjjn. We need to parametrized it as arising from three contributions: 

^A.^ (Af7^)2 + (Ap,)2 ^1 (M7^)2 + (ApT)2/+(AfJ^)2 + (Aft)^' PI (A/J^)2 + (Apt)M' ^ ' 

where the coefficients are the integrated numbers of QED2 photons of mass AI^j, produced per event. 

The small transverse momentum uncertainties and the extension of the experimental data down to small values of 
Pt in pp collisions make it useful to examine first the transverse momentum distribution in pp collisions. In Fig. 1 we 
show the experimental normalized dN^/dpT of anomalous soft photons per pp collision event at 450 GeV/c from Fig. 
2b of Ref. p, [fi^l, after subtracting the bremsstrahlung contributions. The dN'^ /dpr data for pp collisions in Fig. 1 
can be explained by assuming three anomalous soft photon contributions with parameters 

A'/Jy = 17 MeV, = 45 MeV, AfJ^ = 5 MeV, = 0.040, N'{ = 0.017, and = 0.021. (54) 

We examine next the anomalous soft photon transverse momentum distribution in e+-e~ annihilations Q- With an 
error of ApT as large as 10 MeV, the distribution cannot be sensitive to the component at M^t ~ 5 MeV. We show in 
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FIG. 2: (Color Online) Solid points give the DELPHI experimental anomalous soft photon dN^ /dp^ data in hadronic decay 
in e^-e~ annihilations, after subtracting the bremsstrahlung contributions Qj- The solid curve is the total theoretical dN'' /dp^ 
with contributions from component at 13 and 55 MeV, whose separate contributions are shown as the dashed and dashed-dot 
curves, respectively. The dashed-dot-dot curve gives the distribution for the lower momentum component at A/Jj,=5 MeV, 
which however cannot be resolved in the present DELPHI measurements. 

Fig. 2 the DELPHI experimental dN''/dp'^ data from Fig. 4f of Ref. 0, after subtracting the inner bremsstrahlung 
contributions. The experimental data can be explained by assuming the following set of parameters, 

A/o\, = 13 MeV, = 55 MeV, N^^ = 0.106 and A^7 = 0.058. (55) 

There are no reliable data points at < 0.001 (GeV/c)^ to fix the Mjj, component with the present e~^-e~ data. 
For illustrative purposes, we show the M^j- component calculated with Mjy=5 MeV and /Nq =0.525 (as in Eq. 
([M])) shown as the dashed-dot-dot curve in Fig. 2, to indicate that its presence or absence has little effects on the 
theoretical results above pt > 0.001 GeV^. 

Our comparison of the anomalous soft photon transverse momentum distributions reveals that it is necessary to 
examine the transverse momentum distributions of both pp collisions and e~^-e~ annihilations as complimentary 
data sets, as the pp data have finer resolution and smaller errors in the lower pT^15 MeV regions while the e~^-e~ 
data have less fluctuations in the higher pT^50 MeV region. The combined analysis indicates that the transverse 
momentum spectrum can be qualitative described by two components with transverse masses of ~15 MeV and ~50 
MeV, in approximate agreement with the gross features of the theoretical QED2 photon model. There is however an 
additional, lower momentum component at ~5 MeV which shows up in pp collisions, but cannot be resolved in e~^-e~ 
annihilations. The origin of this \ow-pT component is not known and will be left for future studies. Among many 
possibilities, it may be the manifestation of the zero mode of QED2 photon production whose investigation will be of 
great future interest. 

VII. RATES OF MESON AND ANOMALOUS SOFT PHOTON PRODUCTION 

We shall now examine the remaining feature concerning the rates of meson and anomalous soft photon production 
in high-energy e+-e~ annihilations to complete our comparison of the QED2 model with experimental data. There 
are many important physical quantities in the production processes. The receding quark and antiquark generate a 
QCD field of strength k'^^ and a QED field of strength k^^ between the quark and the antiquark in the flux tube 
that will produce the QCD mesons and the QED2 photons, respectively. Here, we have used the superscript a=h 
for hadron quantities and a="f for photon quantities. Each of the field quanta is produced in a final state possessing 
a transverse momentum, and thus the mass that enters into the consideration of quanta production should be the 
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transverse mass Mfj, = \J [Mf^^^^Y +Pt' '^hich can be identified as the boson mass Af"j2D) two-dimensional 

space-time as discussed in Section V. 

To obtain an estimate, we can rely on the Schwinger mechanism of particle production in a strong field as a 
guide (63l - [66j . The probability of particle production of a composite particle of transverse mass Mfrp depends on the 
exponential factor of exp{— 7r(M"j./2)^/K;^g}, where the factor of 1/2 in M"y/2 is to denote the production of a pair 
of particles each of which has a mass Mfrp/2, and the binding of one particle of mass M"j,/2 with the neighboring 
particle of mass M"j,/2 leads subsequently to a composite stable boson of mass Mjrp/2 + Mfj,/2. Furthermore, from 
dimensional analysis, we can infer that the rate of production per space-time volume element {dz dt) has the dimension 
K^^. We therefore assume that the rate of the production of the number of particle of type a, isospin /, and mass 
Mfrp due to the presence of the QCD and QED fields between a receding quark q and an antiquark q is 

^=^E^..-<.-p|-^^|, « = 7,^, (56) 

where Pqq is the probability for the quark-antiquark source pair to be a uu or dd pair, and A is a dimensionless 
constant. In an e^-e" annihilation at high energies, there is an equal probability for the quark-antiquark pair to be 
a uu or dd pair, and so Puu = Pdd — 1/2- 

For the production of QCD2 mesons, the color electric field strength between the leading quark and antiquark is 
independent of the quark flavor quantum number. It is given by 

^tu - = .9QCD2/2 = b. (57) 

For the production of QED2 photons, the electric field strength between the leading quark q and antiquark q is given 
in terms of the electric charges of the quark and the antiquark as 



K 



7 _ in n-u2 
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IQ9QdeQED2/2. (58) 
Thus, wc find that between a receding q and q, there is a constant electric field with a strength 

= 0.002048 GeV^ (59a) 
= 0.000512 GeV^ (59b) 

The QCD field strength and the experimental meson transverse masses as given in Eqs. ((5T|) and ((52|) allow us 
to determine from (|56|) the number of mesons (in a particular Iz state) produced in a space-time volume of AzAi. 
Similarly, the QED field strength k^^, the QED2 photon isoscalar photon mass of M'^j.^Vb MeV, and the isovector 
photon mass of A/^j,~^50 MeV from Eqs. ([5l)) and ([55]) allow us to determine the number of photons produced. We 
obtain, 

= AAzAt X 0.05368 GeV^ (60a) 
= Al\zbd X 0.01391 GeV^ (60b) 
7V7 = AAzAt X 0.0003980 GeV^ (60c) 
Nl ^ AI\zM X 0.0011206 GeV^. (60d) 

In these estimates, the number of produced particles of different types and isospin quantum numbers are proportional 
to the same space-time volume AzAt. This space-time volume fluctuates in each hadronic decay event; the number 
of mesons and photons of different isospin quantum numbers will vary from event to event. However, because all these 
particles in each event are produced simultaneously by the fragmentation of the same q-q string and the same space- 
time volume, the ratio of different isospin spin quantum numbers and types of particles can therefore be proportional, 
on an event-by-event basis. The results in Eq. (|60p give 

4 ' 4' ^ ' 

which reveal that isoscalar photons are more preferentially produced than isovector photons whereas isoscalar mesons 
are much less likely produced than isovector meson (in a particular I3 state). The theoretical ratio of /N^=2.8 
compares approximately well with the experimental ratio of Nq /Nl^2A and 1.8 in Eq. ([51|) and ([55]) extracted by 
fitting the experimental transverse momentum distribution data of [gj and [jj • Using these results, we can also 
construct the ratio of ratios, 
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which states that the number of soft isoscalar photons associated with the isoscalar meson production are more 
numerous than soft isovector photons associated with the isovector meson production. 

The DELPHI experimental measurements Q provide information on the ratios of anomalous soft photon pro- 
duction with various produced neutral or charged meson multiplicities. To compare with experimental data, we need 
to convert the number of different species of mesons to the number of charged and neutral mesons. Isovector mesons 
are pions which have two charged states and one neutral states, and each isoscalar 77 meson decays into 1.64 neutral 
particles (with 2 7's counted as a tt'^ as in Ref. Q) and 0.57 charged particles. Thus, the total meson particle number 
is A'par = (1.64 + 0.57)Nq + 3N^. From Eq. (|60p . the theoretical ratio of total soft photons to total meson particles 
(charged and neutral) is 



7.91 X 10~^ (63) 



which compares reasonably well with the experimental ratio of A^'^'/A^par ~ 9.1 x 10""^. 

In our QED2 photon model, production of mesons of isospin quantum number / will be associated with the 
production of QED2 photons of the same isospin quantum number /. Thus, isoscalar QED2 photons will be associated 
with isoscalar mesons while isovector QED2 photons will be associated with isovector mesons. 

From Eqs. (|60d[) and (|60bp . the theoretical ratio of the number of produced isoscalar photon to the number of 
produced isoscalar meson is 

= 80.6 X 10"^ (64) 

As each isoscalar meson produces 1.641 neutral 7r°-like particles and 0.57 charged particles, the isoscalar meson is 
associated with the production of dominantly neutral mesons. As we consider the production of isoscalar mesons to 
be associated only with the production of isoscalar photons, then for the isoscalar mode of production, N'' /N^cu ^ 
iVp/iVneu which leads to A^^/iVneu ~ -^0 /(1-641A^(5') after summing over all charged particles. From Eq. (|64| . the 
theoretical ratio of soft photon to neutral particle number is 

— = 49.1 X 10~^ 65 

A^neu 1.641A^^' ^ ' 

which comes close to the experimental ratio of N'' /N^cu ^ 37.7 x lO^'^. 



TABLE II: Comparison of QED2 photon model description of the anomalous soft photon production with quantities measured 
or extracted from the pp collision [y| and DELPHI e"'"-e~ annihilation experimental data [3-[3l 



Quantities 


QED2 Model 


Experimental Anomalous Soft Photon Data 


Isoscalar photon mass Mo 


0(25 MeV) 


~15 MeV 


Isosvector photon mass Mi 


0(44 McV) 


~50 MeV 




11/4 


1.8-2.4 




7.91 xlO"-" 


9.1x10"'' 


NyNncu 


49.1x10"^ 


37.7x10"^ 




3.71x10"'' 


6.9x10"'' 



From Eqs. (|60cp and (|60ap . the theoretical ratio of the number of produced isovector photon to the number of 
produced isoscalar meson in a particular state is 

—\- ^ 7.41 X 10"^ (66) 

The isovector meson is three-fold degenerate with two charged particles and one neutral particle. Thus, the production 
of an isovector meson is associated with the production of dominantly charged particles. In the QED2 photon model, 
the sources of isospin current disturbances that produce the mesons and photons are the same. Therefore, the 
production of isovector mesons is associated only with the production of isovector photons. Consequently, we have 
N'^ /Nch N^/Nch, which leads to N'^ /Nch N7/{2Ni) after summing over all neutral particles. Equation ([66)) 
then leads to 
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which is shghtly less than the experimental ratio N^/Nch ^ 6.9 x 10"'^, but is within the same order of magnitude. 

As a summary, we give the comparisons of various quantities obtained in the QED2 model with the DELPHI 
data in Table II. We conclude from the comparison that gross features of the DELPHI data is approximately 
consistent with the QED2 photon model. 



VIII. FURTHER EXPERIMENTAL TESTS OF THE QED2 PHOTON MODEL 

While the QED2 photon model appears to explain qualitatively the main features of the experimental anomalous 
soft photon data, it is desirable to carry out further experimental measurements to test the model: 

1. It will be of interest to measure the transverse momentum distribution of the soft photons with a finer px 
resolution and greater precision for a given narrow range of photon rapidities. Qualitatively, we expect the 
production of photons with two different average transverse momenta, one at ~15 MeV for the production of 
the isoscalar photon and one at ~50 MeV for the production of the isovector photon. 

2. It will be necessary to confirm the presence of the low momentum component at A/Jj,=5 MeV in high-energy 
e+-e^ annihilation experiments. As the origin and the properties of such a low pT component is still unknown, 
additional experimental information on this source of anomalous soft photons will improve our understanding 
of such a component. 

3. It will be of interest to measure the transverse momentum distribution by selecting events with predominantly 
neutral particles and events with predominantly charged particles. The former events will likely arise from the 
production of isoscalar mesons and QED2 isoscalar photons, with an average photon transverse momentum of 
^15 MeV, while the latter from the production of isovector mesons and the isovector photons, with an average 
photon transverse momentum of ~50 MeV. 

4. The rapidity distribution of the produced photons should exhibit the plateau structure, as expected of similar 
distributions in meson production. A measurement of the rapidity distribution will provide useful additional 
information on the dynamics of soft photon production. 

5. Measurements of the properties of associated hadrons similar to those of the DELPHI Collaboration should be 
carried out with hadron-hadron collisions at high energies where anomalous soft photon production has been 
reported [IHll. 



IX. CONCLUSIONS AND DISCUSSIONS 



A color flux tube is formed when a quark and an antiquark (or a diquark) pull apart from each other at high 
energies. The motion of the quarks in the underlying vacuum of the flux tube generates color charge oscillations 
which lead to the production of mesons. As a quark carries both a color charge and an electric charge, the color 
charge oscillations of the quarks in the vacuum are accompanied by electric charge oscillations, which will in turn lead 
to the simultaneous production of soft photons during the meson production process. 

To study these density oscillations, we start with quarks interacting with both QCD and QED interactions in 
four-dimensional space-time in the U(3) group which breaks into the color SU(3) and the QED U(l) subgroups. 
Specializing to particle production at high energies, we find that the dominance of the longitudinal motion and trans- 
verse confinement lead to the compactification from QED4xQED4 in four-dimensional space-time to QCD2xQED2 
in two-dimensional space-time, with the formation of the fiux tube. In the flux tube, the self-consistent coupling of 
quarks and gauge fields lead to color charge and electric charge oscillations that give rise to stable QCD2 bosons and 
QED2 bosons. The boson masses depend on the gauge field coupling constants. The presence of the fiavor degrees of 
freedom leads to isospin dependence of the boson masses, with the isovector meson mass smaller than the isoscalar 
meson mass, but the mass ordering is reversed for the isoscalar photon and the isovector photon. 

As QCD2 and QED2 bosons are stable in the fiux tube environment, we can infer from the quantum field the- 



ory description of particle production in Ref. [21[ that these QCD2 mesons and QED2 photons will be produced 
simultaneously in q-q string fragmentation. Under the condition of adiabaticity with no change of the particle en- 
ergy and longitudinal momentum after the produced particle emerges from the production region, the boson mass in 
two-dimensional space-time turns into the boson transverse mass in four-dimensional space-time. 

The QED2 photon model can explain various features of the anomalous soft photon phenomenon. Because both 
color charge oscillations and electric charge oscillations arise from the same density oscillations of the quarks in 
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the vacuum, both QCD2 meson and QED2 photon will be simultaneously produced by the frag mentation of the 
q-q string. These features are in agreement with those observed in DELPHI experiments [8l 4ld |. The transverse 
momentum distributions of anomalous soft photons in pp collisions @ and e'^-e" annihilations [3 can be described 
by a component with Mqj,^15 MeV and a component at Af^j,~50 MeV in approximate agreement with theoretical 
estimates of the order of the isoscalar and isovector QED2 photon masses. 

In the QED2 model, there are important and non-trivial isospin dependencies in the rate of photon and hadron 
productions that is consistent with recent DELPHI data. The model predicts that the isoscalar photon mass is lower 
than the isovector photon mass. Consequently, the production of isoscalar photons is more likely than isovector 
photons. In contrast, the QCD isoscalar meson mass is greater than the isovector mass, the production of isoscalar 
mesons is less likely than isovector mesons. Thus, the ratio of Nq /Nq can be much greater than Nj/N^. The 
production of isoscalar hadrons is associated with the production of isoscalar photons and leads predominately to 
neutral particles while the production of isovector hadrons is associated with the production of isovector photons 
and leads predominantly to charged particles. As a consequence, the ratio N'^ /N^eu is much greater than the ratio 
TV/A^ch, as observed by the DELPHI Collaboration 

Although the QED2 photon model appears to be promising, it is desirable to carry out additional experimental 
measurements to test the model. We suggest the search for the two components of transverse momentum distributions 
by making appropriate cuts in soft photon rapidities and selecting different regions of neutral and charge multiplicities 
where different isospin photon components are expected. The identification of the two components of different soft 
photon transverse momenta will be a crucial test of the QED2 photon model in the QCD string fragmentationprocess. 

Our examination of the transverse momentum distribution of anomalous soft photons in pp collisions in [6[ reveal 
the presence of an additional component characterized by a transverse mass of M2j,=b MeV. What is the nature 
of this component? Is it related to the zero mode of density oscillations? How does the zero mode manifest itself 
experimentally? Experimental investigation of the low transverse momentum component of photon production and 
the theoretical investigation of the zero mode of QED2 photon production will be of great future interest. 

There are puzzling elements of the QED2 photon model that call for future theoretical and experimental resolution. 
As it now stands, the theoretical determination of the QED2 photon masses is rather uncertain as the mass scale in 
the bosonization of the scalar density is unknown. The QED2 photon mass scale as extracted from experimental data 
requires an electromagnetic current quark mass smaller than the current quark mass as determined from perturbative 
QCD. Whether or not such a smaller value of the mass scale /i for the anomalous soft photon production is justified 
will require further theoretical and experimental investigations. 

Another puzzling and unresolved question is the more detail description of the evolution from a QED2 photon to 
a QED4 photon. We have used the concept of adiabaticity so that the photon preserves its energy and longitudinal 
momentum, only to develop a transverse momentum to balance the mass shell condition. Such a description appears 
to give a qualitative description of the transverse momenta and production probabilities of the soft photons. However, 
there is no additional content in the dynamics of the evolution in our hypothesis. A more detail dynamics of the 
evolution of QED2 to QED4 will be of great interest. 

Finally, if the model is proved to be successful in explaining the anomalous soft photon data, it may be useful to 
explore whether one can study this non-perturbative problem in the full four-dimensional space-time without resorting 
to the intermediate stage of going through the two-dimensional space-time, where non-perturbative physics can be 
carried out more readily. The success of a completely four-dimensional description will provide new insight into the 
non-perturbative behavior of particle production in strong fields. 
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